Stability of two-dimensional ion-acoustic wave packets in quantum plasmas 
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The nonlinear propagation of two-dimensional (2D) quantum ion-acoustic waves (QIAWs) is studied in a 
quantum electron-ion plasma. By using a 2D quantum hydrodynamic model and the method of multiple scales, 
a new set of coupled nonlinear partial differential equations is derived which governs the slow modulation of 
the 2D QIAW packets. The oblique modulational instability (MI) is then studied by means of a corresponding 
nonlinear Schrodinger equation derived from the coupled nonlinear partial differential equations. It is shown that 
the quantum parameter H oc h, associated with the Bohm potential, shifts the MI domains around the fc0-plane, 
where k is the carrier wave number and is the angle of modulation. In particular, the ion-acoustic wave (IAW), 
previously known to be stable under parallel modulation in classical plasmas, is shown to be unstable in quantum 
plasmas. The growth rate of the MI is found to be quenched by the obliqueness of modulation. The modulation 
of 2D QIAW packets along k is shown to be described by a set of Davey-Stewartson-like equations. The latter 
can be studied for the 2D wave collapse in dense plasmas. The predicted results, which could be important 
to look for stable wave propagation in laboratory experiments as well as in dense astrophysical plasmas, thus 
generalize the theory of MI of IAW propagations both in classical and quantum electron-ion plasmas. 

PACS numbers: 52.27.Aj; 52.30.Ex; 52.35.Fp; 52.35.Mw; 52.35.Sb 



I. INTRODUCTION 

The importance of quantum effects has been recognized 
over the last few years in view of its remarkable applications 
in metallic and semiconductor nanostructures (e.g., metal- 
lic nanoparticles, metal clusters, thin metal films, nanotubes, 
quantum well and quantum dots, nano-plasmonic devices etc.) 
OlBE&BfilZllil] as we ^ as m dense astrophysical envi- 
ronments (e.g., white dwarfs, neutron stars, supernovae etc.) 

It is well known that the nonlinear propagation of wave 
packets in a dispersive plasma medium is generically sub- 
ject to amplitude modulations due to the carrier wave self- 
interaction, i.e., a slow variation of the wave packet's enve- 
lope due to nonlinearities. Under certain conditions, the sys- 
tem's evolution may thus undergo a modulational instability 
(MI), leading to energy localization via the formation of en- 
velope solitons. Such solitons are governed by a nonlinear 
Schrodinger (NLS) equation where the nonlinearity at the first 
stage of the amplitude is balanced by the group dispersion. 
This mechanism is encountered in various physical contexts, 
including pulse formation in nonlinear optics, in material sci- 
ence, as well as in plasma physics. A number of works can 
be found in the literature for the investigation of MI in clas- 
sical (see, e.g., Refs. HI H H [3 [ll [MES) as well as 
quantum plasmas (see e.g., Refs. Hm Jl9| lioTl^lll). The MI 
of ion-acoustic waves (IAWs) has been shown to be a general 
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property in a nonlinear dispersive plasma medium, when the 
modulation is considered obliquely to the direction of propa- 
gation of the wave vector [11]. Experimental observations of 
the MI of IAWs have been reported by Watanabe 11211 . Sta- 
ble wave propag ation from modulational obliqueness, in both 
classical |jl3lll4l[l6l[T7ll and quantum plasmas f20ll . has also 
been investigated by a number of authors. 

On the other hand, in a wide variety of scientific fields (e.g., 
nonlinear optics, plasma physics, fluid dynamics etc.), certain 
nonlinear governing equations often exhibit important phe- 
nomena other than solitons, such as shocks (wave singular- 
ity), self-similar structures, wave collapse (i.e., blow up with 
growing amplitude or amplitude decay at finite time or finite 
distance of propagation), as well as the wave radiation emis- 
sion leading to the onset of chaos. One such important system 
in this context is the Davey-Stewartson (DS)-like equations 
[221, where the system has both quadratic as well as cubic 
nonlinearities. Such equations can appear not only in the field 
of fluid dynamics as in the case of water wave propagation 
I22I l23tl . or in optical communications and information pro- 
cessing in nonlinear optical media J24l I25I l26ll . but also in 
plasma physics community I27I l28ll . The DS description of 
collective excitations in Bose-Einstein condenstates has also 
been studied in the context of matter-wave solitons H29l I30TI . 

The purpose of the present work is to investigate, in more 
detail, the stability and instability criteria for the modula- 
tion of quantum IAW (QIAW) packets using two-dimensional 
(2D) quantum fluid model, and to generalize the previous in- 
vestigation I 111], both from classical and quantum points of 
view. Here we show that the nonlinear dynamics of QIAWs as 
well as the static zeroth harmonic field is governed by a new 
set of coupled nonlinear partial differential equations, which, 
in particular, reduces to a set of Davey-Stewartson (DS)-like 
equations for surface water waves [22]. We show that the 
quantum parameter H fr, associated with the Bohm poten- 
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tial, shifts the MI domains in the kO -plane in 2D quantum 
plasmas. Here k is the carrier wave number and 9 is the an- 
gle of modulation with the propagation vector. We also show 
that the maximum growth rate of MI can be reduced by the 
obliquenes parameter 9 rather than H as in one-dimensional 
(ID) quantum plasmas lfl8l[T9tl . Moreover, a criterion for the 
existence of 2D wave collapse is presented from the DS-like 
equations. 

The paper is organized as follows. In Sec. II, we describe 
the 2D quantum hydrodynamic model, and derive the govern- 
ing system of equations using the method of multiple scales. 
In Sec. Ill the MI of the QIAWs is studied and the condition 
for the existence of 2D wave collapse is derived. Finally, the 
Sec. IV is left for concluding the results. 



H. FLUID MODEL AND DERIVATION OF THE 
EVOLUTION EQUATIONS 

We consider the nonlinear propagation of QIAWs in a 2D 
quantum electron-ion plasma. The basic normalized set of 
equations in a two-component unmagnetized quantum plasma 
reads 



^ + v.(»,v) = o, 



dt 
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where V = {d/dx, d/dy) , n e i{\ is the electron (ion) number 
density normalized by their equilibrium value no, V = (u,v) 
is the ion ve locity normalized by the ion-acoustic speed c s = 
y/ 'ksTFe i 'mi with ks denoting the Boltzmann constant, m, the 
ion mass, 7> e = h 2 (37t 2 no) 2 ^ /2kB>n e the electron Fermi tem- 
perature, and n the scaled Planck's constant. Also, H = 
HcOpe/ksTFe is denoting the ratio of the 'plasmon energy den- 



sity' to the Fermi thermal energy, where co P j — y/n^e 1 j&^mj 
is the plasma frequency for the y'-th particle. Moreover, (j) is 
the electrostatic potential normalized by ksTp e /e.The space 
and time variables are respectively normalized by c s / co P i and 
the inverse of CO p i. The electron pressure gradient (V/? e ) and 
quantum force in Eq. (O appear due to the electron degen- 
eracy in a dense plasma with the Fermi distribution function. 
The former can be given (since the equilibrium pressure is 
truly three-dimensional) by the following equation of state 
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where Vp e = \J ksTfe / m e is the Fermi thermal speed of elec- 
trons. 

In order to obtain the evloution equations for QIAW pack- 
ets, we employ the standard multiple-scale technique (MST) 
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in which the space and the time variables are stretched 
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where e is a small parameter representing the strength of the 
wave amplitude and v g = (v gx ,v g y) is the normalized (by c s ) 
group velocity, to be determined later by the compatibility 
condition. The dynamical variables are expanded as 

nj = l + Ze" £ nf{^, r], T)exp[i(k -v- cot)l], (T) 
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and (j)} satisfy the reality condition 



S j = S)" 1 * with asterisk denoting the complex conjugate 
of the corresponding quantity. Notice that in Eqs. ©-(O 
the perturbed states depend on the fast scales via the phase 
(k-r— COt) (where k and co respectively denote the wave vec- 
tor and wave frequency), whereas the slow scales only enter 
the l-th harmonic amplitude. We suppose that k makes an 
angle 9 with the x-axis, and the modulation is along any di- 
rection in the xy-plane. In a general way, the wave vector k 
is then k = (k x ,k y ) = (A: cos 9,ksin 9) and the group velocity, 
Vg = (v gx ,Vgy) = (v g cos9,VgSm9). 

Substituting the expressions from Eqs. ©-(O into the Eqs. 
(Q]i-(|4]i and collecting the terms in different powers of £ we 
obtain for n = 1,1 = 1 the linear dispersion relation for the 
normalized wave frequency (co — > Co/cOpi) and the wave num- 
ber (k — > kcJcOpi) as 



CO 



k 2 A 
l+k 2 A' 
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where A = 1/3 +H k /4. Note that the wave number k in 
Eq. ( [Tot is not very small such that the wave frequency co 
is sufficiently large to prevent the appearance of harmonic 
modes of k x ,k y as proper modes. These harmonic modes 
will virtually appear in higher orders of perturbations. Also, 
since k is normalized by the inverse of Fermi screening length 
{Xf = c s / COpi), the values of k greater than unity is inadmis- 
sible, otherwise the wavelength would become smaller than 
the screening length. As a result, the collective behaviors of 
the plasma will disappear. Moreover, we consider the quan- 
tum parameter H, to vary in the range 0.1 < H < 0.45 such 
that the ratio Vp>/c <C 1 [an approximate condition for the 
nonrelativistic quantum hydrodynamic model to be valid] and 
the coupling parameter, gQ = 2m e e 2 / 'eon 2 (37Z 2 \fno) 2 l^ < 1 
(which corresponds to the density region where the quantum 
collective and mean field effects become important). 

On the other hand, for the second order reduced equations 
with n = 2,1 = 1, we obtain the following compatibility con- 
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ditions for the group velocity components 



in 2D quantum plasmas 
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Next, proceeding in the same way as in Refs. HH [l9l l2lll 
and finally considering the equations for / = 1 and n = 3, we 
obtain the following coupled equations for the propagation of 
modulated QIAWs 
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>V = *o 

Jd n v^dp X = Jd^\ 2 d7],X = Id n \<l>\ 2 dp and the 
coefficients Pj ,Qft,Rj,Sfi are given in Appendix A, where 
j = 1,...,6 for Q's, j — 1,2,3 for P's and j — 1,...,4 for 
others. The corresponding coefficients Pj ,QpR C j,S c j for 
the propagation of 2D classical IAWs can also be obtained 
by the similar method as discussed in Appendix B. The 
superscripts 'q' and 'c' are used to denote the coefficients 
corresponding to 2D quantum and 2D classical electron-ion 
plasmas. Thus, we have obtained a new system of nonlocal 
nonlinear equations, which describe the slow (and general) 
modulation of the QIAW packets in 2D quantum plasmas. 
The coefficients Pf ,P| appear due to the wave group dis- 
persion and 2D motion, and P% for the arbitrary orientations 
of the carrier wave propagation as well as the modulation of 
the QIAW packets. The nonlinear coefficients Qf is due to 
the carrier wave self-interaction originating from the zeroth 
harmonic modes (or slow modes), i.e., the ponderomotive 
force, and Q q 5 ,Ql are for the combined effects of 2D motion 
and self-interaction. The nonlinear-nonlocal coefficients 
62 >2f >Ql al "i se due to the coupling between the dynamical 
field associated with the first harmonic (with a 'cascaded' 
effect from the second harmonic) and a static field generated 
due to the mean motion (zeroth harmonic) in the plasma. The 
appearance of the coefficients RpS?- can also be explained 
similarly. 

Equations (fT3l and < TT~4T > can be studied, as for example, for 
the modulation of Stokes wave train (plane wave) with con- 
stant amplitude to a small perturbation as well as to look for 
envelope solitons, and wave collapse, if any. There are, of 
course, many other physical insights (e.g., dromion solution), 
that can also be recovered from this system. In particular, 
looking for the modulation of QIAW packets parallel to the 
carrier wave vector, one can obtain the DS-like equations 



+ p ?j£ + p 2j^ + Qu\^ + QhM = o,(i5) 
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where the coefficients are those obtained at 9 = from the 
general expressions. Also, Q q n = Q\ + QiSf/Ri and £?22 — 
62 ^ Qi R i/ R t Similar forms like Eqs. O and fl6) can also 
be obtained by transforming the x-axis through an angle 9 , so 
that the wave vector k and the modulation direction coincide. 
Now, our aim is to investigate the MI of QIAWs by consid- 
ering its modulation along the x-axis through an angle 9 with 
the carrier wave vector k. To this end, we disregard the group 
velocity component of the modulated wave along the y-axis, 
as well as the y- or rj- dependence of the physical variables. 
Thus, we obtain from Eqs. ( fl3l ) and ( TT4l > the NLS equation for 
the oblique modulation of QIAW packets 



dT dP 



Q\<p\ 2 $ = 0, 



(17) 



where P = P c ( and Q = Q\ + Q\S\lR\ 



III. ANALYSIS OF MODULATION INSTABILITY 

We consider the MI of a plane wave solution of Eq. ( fT7b 
for (j) with constant amplitude <pQ. The boundary conditions 
that — > as S, — > °o must now be relaxed, because the plane 
wave train is still unmodulated and the solution is not unique. 
Thus, we can represent the solution as the monochromatic 
solution = </>oexp(/g|0o| 2 T), where A(t) = — 2l0o| 2 is the 
nonlinear frequency shift. To study the stability of this so- 
lution we modulate the amplitude against linear perturbation 
as (j) = [0o + m cos(^T4 — ilT)]exp(/g|0o| 2 l : ), where K and Q. 
are, respectively, the wave number and the frequency of mod- 
ulation. We then readily obtain from Eq. ( fTTI i the dispersion 
relation 



p2 K 4 



K 2 



(18) 



where K c = \/2(2/P|0o| is th e critical wave number such 
that the MI sets in for a wave number K < K c , i.e. for all 
wavelengths above the threshold, X c = 2n/K c . The instability 
growth rate (letting Q. = iF) is then given by 



t=pk* 



K 2 



1. 



(19) 



The maximum growth rate, r max = |2ll0o| 2 as equals to the 
amount of the nonlinear frequency shift |A(t)|, is achieved at 
K = K c / v2- Clearly, the instability condition depends only 
on the sign of the product PQ, which, in turn, depends on 
the angle of modulation 9 as well as the quantum parameter 
H. These may be studied numerically, relying on the exact 
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expressions of P and Q. We find that for a fixed value of 
H = 0.3 (i.e., as one enters the density region no = 1.5295 x 
10 33 m- 3 ), P > in 0.6958 < < n/2 and Q < in < < 
0.5281 and 1.138 < < 7Z;/2 for < /t < 1. Other domains of 
k and where P and Q change sign depend on the P = and 
Q = curves in the £0-plane. Similarly, considering another 
value of H = 0.45 (where n = 1.34 x 10 32 m~ 3 ) one can find 
P > in 0.6748 < < n/2 and g < in < < 0.6 and 
1.137 < < n/2 for < k < 1. As above, P,g also change 
their sign in other regions in the fc0-plane. 

Figure 1 displays the contour plots of PQ = boundary 
curves against the normalized wave number k and the mod- 
ulation angle 0. The upper panel corresponds to the quan- 
tum case, whereas the lower one is for the classical electron- 
ion plasmas. The stable (PQ < 0) and unstable (PQ > 0) re- 
gions are separated in the kO -plane by the boundary curves of 
PQ = as indicated by the white and the gray regions respec- 
tively. We have allowed k to vary in a range < k < 1 (as 
explained in Sec. II) and to vary between 0=0 and = n, 
so that all plots seem to be n/2 periodic, and, in fact, sym- 
metric upon reflection with respect to either 0=0 or 9 = n/2 
lines. We also consider H, in the range 0.1 < H < 0.45 such 
that gg < 1 (as explained before). It is clear from Fig. 1 that 
the stability and instability domains are quite different in clas- 
sical and quantum plasmas. Evidently, the quantum parameter 
H shifts the MI domains in the £0-plane. It is, of course, hard 
to find the common instability regions for both classical and 
quantum cases. Increasing the value of H (< 0.45) gives no 
effective change, except in reducing the stable regions around 

In particular, considering the modulation of the QIAWs 
along the carrier wave vector k (i.e., the case of = 0), the 
QIAW is shown to be modulational unstable (see Fig. 2). 
This is in contrast to the classical case where the IAWs are 
known to be stable under parallel modulation [11]. Physically, 
the QIAW under parallel modulation becomes unstable due to 
the nonlinear self-interactions originating from the zeroth har- 
monic modes or slow modes as well as the second harmonic 
modes (since P and Q are always negative). However, the case 
of modulation perpendicular to the wave vector k (0 = n/2) 
gives rise stable wave propagation, an agreement with clas- 
sical plasmas 111 111 . Here the wave is stable due to the second 
harmonic self-generation of the modes, since P > and Q < 0. 
Thus, for the MI to occur the second harmonic mode is very 
essential in order to maintain the same sign of P and Q. 

On the other hand, the MI growth rate as given by Eq. ( fT9] i 
can be calulated as depicted in Fig. 3. The upper panel (cor- 
responding to QIAWs) shows that the growth rate can be re- 
duced by increasing the obliqueness of modulation giving rise 
cut-offs at lower wave numbers of modulation. The quantum 
parameter H has no significant role in reducing such growth 
rate. In the classical case, the maximum growth rate reduces 
slightly with a small change of the angle of modulation (see 
lower panel of Fig. 3) 
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FIG. 1: The equation PQ = is contour plotted in the k9 -plane to 
show the stable and unstable regions in quantum (upper panel) and 
classical (lower panel) plasmas fEq.d 17t1 for H = 0.3. The white or 
stable (gray or unstable) areas correspond to the regions in the k6- 
plane where PQ<0 (PQ > 0). 




FIG. 2: The case of parallel modulation (6 = 0): PQ is plotted 
against k to show that the QIAW is unstable (PQ > 0). The solid 
and dashed lines respectively correspond to the cases where H = 0.3 
and H = 0.45. 



IV. VARIATIONAL PRINCIPLE AND 2D-EVOLUTION 

The DS-like equations $15[ and ( fTSI l describes the case of 
general 2D-evolution of slowly varying weakly nonlinear ion- 
acoustic waves in the quantum regime. Here we start by point- 
ing out a variational principle for the DS-like system. Intro- 
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principle makes it possible to use trial functions as a means to 
obtain approximate solutions. A good review on this topic for 
wave equations similar to (fl~5T > and ( fToT ) is given by Ref. l39ll . 
A thorough variational study of ( fl"5l > and ( [ToT l for the general 
case is beyond the scope of the present paper. However, in or- 
der to illustrate the usefulness of the variational technique, we 
will state a variational result that applies when a cylindrically 
symmetric collapse is possible. For the general case, Eqs. 
( TTBl l and ( TToT l do not admit a cylindrically symmetrical col- 
lapse. However, for geometries and parameter regimes where 
the coefficients fulfill either R 2 < Ri, or (^ 2 S\/R\ -C Q q n , 
our system reduces to a 2D NLS equation for which cylindri- 
cally symmetrical solutions are possible. Assuming that at 
least one of the above strong inequalities apply, and that we 
have cylindrically symmetric initial conditions (to be defined 
below), the evolution equation for our system can be written 



dtp 



i a 

r dr 



dtp 



(24) 



FIG. 3: The MI growth rate T is shown with respect to the wave 
number of modulation, K. The growth rate can be quenched by the 
obliqueness parameter 6 in quantum plasmas (upper panel). In the 
classical case (lower panel), the maximum growth rate is reduced 
slightly with a small change of 9. The solid, dashed and dotted lines 
correspond to the values 6 = 0.99, 1.0, 1.05 respectively. 



ducing the Lagrangian density 
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where the action functional is £/(<p,tp*,u) = / ^fd£,dz, and 
d 2 u/dt, 2 = y/ plays the role of a potential. We obtain eqs. 
(Q3) and ([ToT l by varying <j), <j>* and u and minimizing the ac- 
tion as usual. It is straightforward to show that $15[ and ( [ToT l 
posses constants of motion, representing the conservations of 
the number of high frequency quanta 



where we have introduced the radial coordinate r 2 = t, 2 ■ JP\ + 
TJ 2 /Pj [40], and assumed, cylindrical symmetry, i.e. that the 
spatial dependence is only through this variable . The coef- 
ficient Qeff is either Q q n (in case the reduction from (fl5l l and 
( fT6] > was due to the inequality 622^1/^1 ^ 2?i) 01 2eff = 
2li + $riS\lR\ ( m case me reduction was due to R2 -C R\)- 
For the case of Q e ff > 0, collapse is possible. A good es- 
timate for the collapse threshold can be obtained by using 
the variational formulation, with a trial function of the form 



0(T,r) = F(T)sech(r//(T))exp(/fr(T)r 2 ), see e.g. Ref. |4JJ] 
for mathematical details. The result is that there will be a 
collapse towards a zero radius of the pulse profile, /(t) — » 
within a finite time, provided the initial profile is sufficiently 
intense and well localized, as expressed by the collapse con- 
dition 



2eff/ 2 (T = 0)|F 2 (T = 0)1 > 1.35 



(25) 



The accuracy of this collapse condition is well supported by 
numerical calculations. Naturally, for cases where cylindrical 
symmetry does not apply, the evolution will be much more 
complicated than indicated by our simple example. 



d 

dx J-o 



\tp\ 2 drid% =0 



(21) 



and of longitudinal momentum 



d_ 

dx 



dr]d% = (22) 



as well as transverse momentum 



— f 

dx J-, 



drjdZ =0. 
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Furthermore, since no explicit time dependence occurs in the 
Lagrangian, the Hamiltonian J™ a ,J™ a ,Jf?drid£, derived from 
Jz? in ( l20l i. is also conserved. The existence of a Variational 



V. CONCLUSION 

We have investigated the nonlinear propagation of QIAW 
packets in a 2D quantum plasma. A multiple scale technique 
is used to derive a coupled set of nonlinear partial differential 
equations, which governs the dynamics of modulated QIAW 
packets. The set of equations, in particular (i.e., modulation 
along the direction of the pump carrier wave), is shown to be 
reducible to a well-known Davey-Stewartson (DS)-like equa- 
tions B22I1 . The latter can be studied for the 2D wave collapse 
in dense plasmas. The oblique MI of QIAWs is then stud- 
ied by means of a corresponding NLS equation. It is found 
that the quantum parameter H has the significant role in shift- 
ing the instability domains around the £0-plane. The case of 
parallel modulation, which is known to be stable in classical 
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plasmas fiUl . is shown to be unstable in quantum plasmas. In 
contrast to the classical case, the MI growth rate in quantum 
plasmas can significantly be reduced by increasing the angle 
of modulation. 

It is to be mentioned that the numerical simulation of the 
DS equations could furnish more convincing evidence for the 
process of soliton collapse, where the quantum parameter H 
may play a crucial role in accelerating or decelarating the col- 
lapse process. However, it needs extra effort, and could be an 
open issue to be explored in future studies. The integrability 
as well as the dromion solution, if there be any, of the system 



will also be an another investigation, but beyond the scope of 
the present work. 
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APPENDIX A 

The coefficients Pj,gJ,i?^,5y appearing in Eqs. (fT3T i and 
(fl4] > are given as follows. 
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In the expression for Q\ the first (second) term in the square 
brackets is due to the first and second harmonic modes (the 
first and zeroth harmonic modes). 
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APPENDIX B 



The coefficients corresponding to classical electron-ion 
plasmas can be obtained by considering a 2D fluid model 
with contininuity and momentum equations for cold ions and 
Boltzmann distributed electrons [see, e.g., the Eqs. (5)-(8) in 
Ref. ifTTIl l. The normalizations and the corresponding fluid 
equations for the variables can be recovered by simply replac- 
ing the Fermi temperature Tf e by the classical temperature T e 
with considering the electron pressure as p e — ksT e n e instead 
of the Fermi pressure, and disregarding the Bohm potential 
term proportional to h. These coefficients are presented for 
a general interest of the readers to study the dynamics of ion- 
acoustic waves from the evolution equations of the forms ( fT3l , 
(fl4l > or (JT3J, ( [ToT l in 2D classical electron-ion plasmas. 

Thus, the coefficients Pj ,Q C j,R C j,S C j to be appeared in equa- 
tions like ( fT3] l and (fT4l > in the case of classical plasmas can be 
given as follows: 
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where, in the expression for Q c v the first (second) term in the 
square brackets is due to the first and second harmonic modes 
(the first and zeroth harmonic modes). Moreover, 
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where ft) 2 — k 2 /(l + k 2 ) is the normalized classical disper- 
sion relation (ft) — > to/co P i and k — > kc s /(0 P i). We note that the 
difference in the dispersion relation as compared to Eq. (10) 
stems from the use of a classical (i.e. not a Fermi) equation of 
state. 
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